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In  this  paper,  we  present  a  robust  control  of  a  variable  speed  Doubly  Fed 
Induction  Generator  (DFIG)-based  Wind  Energy  Conversion  System 
(WECS),  using  Fractional  Order  Control  (FOC)  to  prevent  system 
deterioration  under  different  critical  conditions  (external  disturbance, 
measurement  noise  and  DFIG  parameters  variation).  In  order  to  extract  the 
maximum  power  from  the  wind,  a  Maximum  Power  Point  Tracking  (MPPT) 
strategy  based  on  rotor  speed  control  is  proposed.  Furthermore,  a  vector 
control  strategy  is  used  for  controlling  active  and  reactive  powers  of  DFIG. 
Additionally,  a  simple  design  method  of  Fractional  Order  Proportional 
Integral  (FOPI)  controller  is  proposed.  Finally,  the  system’s  performance  is 
tested  and  compared  according  to  reference  tracking,  robustness,  disturbance 
rejection  and  noise  minimization. 
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1.  INTRODUCTION 

Nowadays,  a  majority  of  the  wind  energy  conversion  systems  utilize  doubly  fed  induction  generator 
because  of  its  high  efficiency  and  low  cost  owing  to  its  configuration  [1],  [2].  DFIG  have  windings  on  both 
stator  and  rotator  parts,  where  both  of  them  transfer  considerable  power  between  generator  and  grid.  In 
DFIG,  the  converters  process  only  about  20-30%  of  total  generated  power,  and  the  rest  is  fed  to  the  grid 
directly  from  the  stator  [3] -[5]. 

Usually,  DFIG  is  controlled  using  vector  control  strategy,  which  is  either  stator  voltage  oriented,  or 
stator  flux  oriented.  This  scheme  decouples  nonlinear  Multiple  Input  Multiple  Output  (MIMO)  system  of 
DFIG  into  two  linear  Single  Input  Single  Output  (SISO)  subsystems  representing  direct  and  quadrate  rotor 
currents  separately  [6]-[8],  The  active  and  reactive  power  control  of  DFIG  is  attained  by  controlling  these 
two  linear  first  order  SISO  subsystems  with  two  PI  controllers.  The  drawback  of  the  vector  control  strategy  is 
that  the  system  performance  depends  on  the  DFIG  parameters,  especially  the  rotor  resistance  [6],  [7],  Thus, 
the  performance  degrades  when  the  parameters  of  DFIG  used  in  the  control  system  design  are  altered  due  to 
temperature  variation,  saturation,  etc. 

To  solve  this  problem,  we  have  replaced  the  conventional  proportional  integral  controller  with  a 
fractional  order  proportional  integral  one  which  is  a  generalization  of  the  former  [9],  [10].  The  fractional 
order  PI  controller  has  a  potential  to  increase  system  robustness  and  improve  the  system  performance  with 
the  additional  parameter  X,  which  is  the  fractional  order  of  the  integration  action  [11],  [12], 


Journal  homepage:  http://iaescore.com/journals/index.php/IJPEDS 


Int  J  Pow  Elec  &  Dri  Syst 


ISSN:  2088-8694 


□  1073 


In  this  work,  our  contribution  to  DFIG  robust  control  is  the  use  of  a  simple  and  practical  fractional 
order  PI  controller,  which  leads  to  good  dynamic  performance  and  robustness  of  DFIG  in  both  healthy  and 
critical  conditions. 

In  section  2,  we  introduce  a  brief  description  of  the  wind  energy  conversion  system  to  be  studied. 
Then,  the  construction  of  a  mathematical  model  for  each  component  (DFIG  and  wind  turbine)  of  the  wind 
system  is  presented  in  section  3.  Next,  the  control  strategy  is  developed  in  section  4.  The  design  of  fractional 
order  PI  controller  is  established  in  section  5  followed  by  the  simulation  results  and  discussion  in  Section  6. 
Finally,  we  complete  this  work  with  a  conclusion. 


2.  SYSTEM  DESCRIPTION 

The  WECS  to  be  studied  is  presented  in  Figure  1 .  This  system  consists  of  a  wind  turbine  driving  a 
DFIG  through  a  gearbox.  The  rotor  is  connected  through  two  converters,  and  the  stator  of  the  DFIG  is 
directly  connected  to  the  grid. 


Figure  1 .  DFIG-based  WECS  structure 


3.  MODELING  OF  THE  SYSTEM 

3.1.  Wind  turbine  model 

In  WECS,  due  to  the  different  losses,  the  extracted  power  available  on  the  rotor  of  the  wind  turbine 
can  be  expressed  as  [13],  [14]: 

Pt  =  \  Cp(TSR,  /?)p  S  v3  (1) 

where  Pt  is  the  turbine  power,  Cp  is  the  power  coefficient,  p  is  the  air  density,  S  is  the  swept  turbine  area  and 
v  is  the  wind  speed.  The  power  coefficient  is  influenced  by  the  pitch  angle,  [3,  and  the  tip-speed  ratio,  TSR  = 
R-S^±i  where  R  is  the  turbine  radius  and  12 1  is  the  turbine  rotational  speed. 

The  wind  turbine  model  is  shown  in  Figure  2,  where  G  is  the  speed  gain  of  the  gearbox,  Taer  is  the 
aerodynamical  torque,  Tg  is  the  gearbox  output  torque,  J  is  inertial  moment,  /  is  viscous  friction  coefficient 
of  the  rotor,  and  ft  is  the  rotor  speed. 

3.2.  DFIG  model 

The  stator  and  rotor  equations  of  the  DFIG  machine  are  derived  from  Park  reference  frame  rotating 
at  synchronous  speed  [15],  [16],  and  can  be  described  by: 

The  stator  and  rotor  voltage  equations: 

—  Ps^sd  T  —<Psd  —  (2) 

^sq  —  PsGq  T  ~(Psq  T  tOs^sd  (3) 

f/'ri  —  Pr^-rd  T  ~(Prd  ~  (ws  —  tif)(prq  (4) 
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Vrq  Rri-rq  T  dt(Prq  T  (ws  w)<Pr 


(5) 


where  Vsd  and  Vsq  are  the  direct  and  quadrate  components  of  stator  voltage  (V),  Vrd  and  Vrq  are  the  direct  and 
quadrate  components  of  rotor  voltage  (V),  isd  and  isq  are  the  direct  and  quadrate  components  of  stator  current 
(A),  ird  and  irq  are  the  direct  and  quadrate  components  of  rotor  current  (A),  Rs  and  Rr  are  the  stator  and 
rotor  resistances  (£2)  and  cjs  and  o>r  are  the  stator  and  rotor  pulsations  (rad/s). 

The  stator  and  rotor  flux  equations: 


* Psd 

^s^sd  Lm  Wd 

(6) 

^Psq 

~  Lsi-sq  Lmirq 

(7) 

(Prd. 

—  Lrird  +  Lmisd 

(8) 

<Prq 

—  Lrirq  +  LmiSq 

(9) 

where  <psd  and  <psq  are  the  direct  and  quadrate  components  of  stator  flux  (wb),  (prd  and  cprq  are  the  direct  and 
quadrate  components  of  rotor  flux  (wb),  Lm  is  the  mutual  inductance  (H)  and  Ls  and  Lr  are  the  stator  and 
rotor  inductances  (H). 

The  electromagnetic  torque  equation  is  given  by: 


Pem  —  P  Lm(irdisq  Wq^-sd) 


(10) 


where  P  is  the  number  of  pole  pairs. 


4.  CONTROL  STRATEGY 
4.1.  Wind  turbine  control 

In  order  to  extract  maximum  power  from  the  wind,  we  apply  the  Maximum  Power  Point  Tracking 
based  on  the  rotor  speed  control.  The  tip-speed  ratio  is  tuned  to  TSR0Vtimai  over  different  wind  speeds  v,  by 
adapting  the  rotor  speed  fl  to  H  ref  —  (G  v  TSR  optimal')/ R  [17].  The  MPPT  scheme  is  shown  in  Figure  2. 


Figure  2.  MPPT  based  on  rotor  speed  control 


4.2.  DFIG  control 

We  apply  vector  control  strategy  to  the  DFIG,  which  decouples  the  DFIG  model  into  two 
independent  subsystems  of  flux  and  torque,  in  order  to  get  the  performance  similar  to  DC  motor  speed 
control  [18],  Thus,  the  stator  flux  is  chosen  to  be  oriented  to  direct  axis  in  the  park  reference  frame  ((psd  =  <ps 
and  <psq  —  0),  and  by  neglecting  stator  resistance,  Vsd  and  Vsq  become:  Vsd  —  0  and  Vsq  —  Vs  —  cos  <ps. 
Consequently,  the  electromagnetic  torque  (Tem  DFIG )  and  the  active  and  reactive  stator  powers  (Ps,  Qs) 
become: 

Tem=  -P<PsL~rirq  dD 

Ls 
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P  =  -V-^i 


rq 


Qs  =  7rr-Vstird 

w5Ps  Li  5 

Also,  the  rotor  voltages  can  be  expressed  as: 

Vrd  =  Rrird  +  (X  -  Lf)  ^  -  slip  0)s  (Lr  -  &)  l 

Vrq  =  Rrirq  +  (Lr  -  t)^T  +  SilP  ^  {Lr  ~  t)  ird  +  Slip 


1 -rq 


LmYs 

Ls 


(12) 

(13) 


(14) 

(15) 


where  slip  is  the  slip  of  DFIG. 

The  proposed  control  plan  of  DFIG  based  on  two  FOPI  controllers  is  shown  in  Figure  3. 


Figure  3.  The  control  scheme  of  DFIG  with  FOPI  controllers 


FRACTIONAL  ORDER  PI  CONTROLLER  TUNING  PROCEDURE 

In  the  previous  section,  the  general  form  of  the  transfer  function  of  the  plants  was: 

P(s)  =  — 

rs+i 


06) 


where  T  and  K  are  constants. 

In  this  section,  two  different  controllers  are  proposed:  the  integer  order  PI  and  the  Fractional  Order 
PI  (FOPI).  The  transfer  functions  of  the  two  controllers  can  be  described  as  follows: 


CP/(s)  =  K„  (l  +  f) 
CF0P,(s )  =  Kv  (l  +  J) 


(17) 

(18) 


where  Kp,  K t  and  X  are  real  numbers,  A  E  [0,1], 


5.1.  Design  specifications 

In  [19],  [20],  a  tuning  method  for  PI  controller  and  FOPI  controller  is  proposed.  We  choose  the  gain 
crossover  frequency,  coc ,  and  phase  margin,  cpm,  to  be  the  same  for  both  controllers.  For  the  system  stability 
and  robustness,  the  following  constraints  are  considered:  [19],  [20] 
a)  Phase  margin  constraint: 


Arg[G(j(oc)\  =  Arg[C(ja)c)P<ja)c )]  =  zC(y'njc)  +  zP(/'wc)  =  -n  +  cpm  (19) 


Robust  Control  of  Doubly  Fed  Induction  Generator  Using  Fractional  Order  Control  (Abdelaziz  Azza) 


1076  □ 


ISSN:  2088-8694 


where  G(jco)  is  the  open-loop  transfer  function  of  the  system,  C(jco)  is  the  controller  transfer  function 
and  P(jtt>)  is  the  plant  transfer  function, 
b)  Gain  cross-over  frequency  constraint: 


\G(jcoc)\dB  =  |C(/nJc)P(/wc)|dB  =  0 


(20) 


c)  Robustness  to  loop  gain  variations,  which  demands  that  the  phase  Bode  plot  to  be  flat  at  the  gain 
crossover  frequency,  a>c,  i.e.,  the  derivative  of  the  open  loop  phase,  at  the  gain  cross-over  frequency,  to  be 
equal  to  zero: 


d(Arg[G(ja>c)]) 


dcj 


=  0 


(21) 


5.2.  Fractional  order  PI  controller  Tuning 

Based  on  the  design  specifications  introduced  in  the  previous  section,  we  present  the  tuning 
procedure  of  the  Fractional  Order  PI  controller,  for  the  first  order  plant. 

The  open-loop  transfer  function  with  the  FOPI  controller  is. 


G(s)  =  CF0PI(s)P(s )  =  Kp  (l  +  J)  (; 


— ) 
Ts+lJ 


where  K  and  T  are  known  and  Kp ,  Kt  and  A  should  be  designed  in  the  controller  design  process. 
The  FOPI  controller  can  be  expressed  as: 


J FOPI 


w  =  *»(!+  ?)  =  *»(!  +  sSl) 


Since  j  =  e n!2,  thus  jA  =  eAn /2  =  cos  (An /2)  +  j  sin(An/2 ),  which  gives, 


- FOPI 


C s )  =  KP  ( 


1  +■ 


cos(An/2)+ j  sin(ln/2)J 

The  open-loop  phase  at  the  gain  cross-over  frequency  is. 


Arg\G  (ja>c)\  =  —arctan 


Ki(x)ca  sin(An/2) 
1+Ki(i>c*  cos(Xn/Z) 

According  to  the  first  design  constraint  (a). 


Arg\G  (/u>c)|  =  —arctan 


sin(An/2) 


l+Ki<jLic  A  cos(Xn/2 ) 

Hence,  the  relationship  between  and  X  can  be  established  as, 

_  -  tan[arctan(  a> CT  )+  <pm] 

1  sin(  X  n  2  )+io^ *  cos(  X  n  2  )  tan[arctan(  wcT  )+<pm] 

And  the  open-loop  gain  using  at  the  crossover  frequency: 


—  arctan(a)cT ) 


arctan(a>cT )  =  —n  +  cpm 


\G(ja>c)\ 


KKp  I[Kicoca  sin(An/2)]2  +  [l+KiQ)cA  cos(An/2 )]2 


V(7Voc)2  +  l 

Additionally,  according  to  the  second  constraint  (b): 

KKpJ[KiWc*  sin(An/2)]2  +  [1+K cos(An/2)]2 

V(T«C)2  + 1  _ 


(22) 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


Also,  according  to  the  third  constraint  (c): 
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d(Arg[G(ja>c)]')\  _  1  sin(l7r/2) _ T  _  g 

dui  <i>^+2Ki<i>c*'  cos(l7r/2)+Jf?  (7’mc)2  +  1 

Moreover,  another  equation  about  KL  can  be  established  as:  Cu>f2AK2  +  BKt  +  C  —  0 

where,  B  —  2Ca>fA  sin(zbr/2)  —  Xu )f*~1  sm(Xu/2)  and  C  —  T 2+1 

Thus, 


Ki  = 


-B±  Bz-4Czcu~ 


2  C(j)r 


(31) 


The  Equations  (27),  (29)  and  (31)  of  the  three  variables  X,  Kt  and  Kp  can  be  solved  easily  by  different 
approaches. 

For  tuning  the  FOPI  controller  and  classical  PI  controller  for  the  rotor  current  control  loop,  the  same 
cross-over  frequency,  a)c,  is  presumed:  ooc  =  500  rad/s.  Furthermore,  in  order  to  get  a  damping  ratio  of  {  — 
0.707,  the  phase  margin  should  be  the  same  for  both  controllers:  tp)n  =  64°.  Hence,  the  tuned  current 
controllers  give, 

CF0PI(s )  =  0.0763  (l  +  (32) 

CP/(s)  =  0.0792  (l  +  (33) 


6.  RESULTS  AND  DISCUSSION 

In  order  to  verify  and  study  the  efficiencies  of  the  control  strategy,  two  sets  of  simulations  were 
performed.  The  first  set  is  carried  out  on  an  isolated  DFIG  current  control  loop  to  prove  the  effectiveness  of 
the  proposed  controller  design.  Whereas,  the  second  one  is  carried  out  on  DFIG-based  WECS  to  confirm  the 
dynamic  of  the  whole  system.  The  WECS  system’s  parameters  used  during  simulations  are  listed  in  Table  1. 


Table  1.  Doubly  fed  induction  generator  and  wind  turbine  parameters 


Doubly  fed  induction  generator  parameters  Wind  turbine  parameters 


Rated  power,  Pn 

1.5  MW 

Blade  radius,  R 

41  m 

Stator  rated  voltage,  Vs 

398/690  V 

Number  of  blades 

3 

Stator  rated  frequency,  fs 

50  Hz 

Gearbox  ratio,  G 

80 

Stator  inductance,  Ls 

Rotor  inductance,  Lr 
Mutual  inductance,  Lm 
Stator  resistance,  Rs 

Rotor  resistance,  Rr 
Number  of  pair  of  poles,  P 

0.0137  H 

0.0136  H 
0.0135  H 

0.012  n 
0.021  n 

2 

Maximal  power  coefficient,  Cp_max 

0.55  (TSRoptimal  =  8.12  and  p  =  0°  ) 

Case  1:  The  isolated  DFIG  (without  WECS)  current  control  loop 

This  simulation  analyses  and  compares  the  performance  of  the  FOPI  controller  with  the  classical  PI 
controller  of  the  DFIG  current  control  loop.  Under  four  different  specifications,  i.e.  the  set-point  reference, 
the  robustness,  the  disturbance  rejection  and  the  noise  minimization. 

Figure  4  shows  the  unit  step  responses  of  the  DFIG  current  control  loop  with  the  FOPI  controller  as 

shown  in  Figure  4(a),  and  classical  PI  controller  as  shown  in  Figure  4(b).  The  unit  step  responses  are 

11  1 

designed  with  three  different  closed-loop  gain  (If):  — ,2 —  and  0.5 — .It  is  clear  that  the  FOPI  controller 

Ry  Ry  Ry 

exhibits  a  good  reference  tracking  and,  on  contrary  to  the  classical  PI  controller,  its  overshoots  were  not 
affected  by  the  different  gains.  Therefore,  the  controlled  system  using  the  fractional  controller  remains  robust 
over  gain  changes. 
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Figure  4.  Step  current  responses:  (a)  with  FOPI  controller,  (b)  with  classical  PI  controller 


The  unit  step  responses  of  DFIG  current  control  loop  with  the  FOPI  controller  and  the  classical  PI 
controller  are  depicted  in  Figure  5.  A  disturbance  of  0.5  is  applied  on  the  output  of  the  plant  from  t  —  0.2  s 
on,  and  a  measurement  noise  with  Gaussian  distributed  random  signal  with  zero  mean  and  0.1  variance  value 
is  assumed  in  the  output  of  the  plant  from  t  —  0.4  s  on.  By  comparing  the  two  curves  we  notice  that  a  faster 
disturbance-rejection,  and  better  minimization  of  the  measurement  noise  effect  is  achieved  using  the  FOPI 
controller,  indicating  superior  disturbance-rejection  and  measurement  noise  minimization  for  the  controlled 
system  with  the  fractional  controller. 


-ft — 

—  Classical  PI 

—  Fractional  order  PI 

7 

J 

I - 1 - 1 - 1 - 1 - 1 - 

0  0.1  0.2  0.3  0.4  0.5  0.6 

Time  (s) 


Figure  5.  Step  current  responses  using  FOPI  controller  and  classical  PI  controller 


Case  2:  the  DFIG-based  wind  energy  conversion  system 

In  this  case,  we  simulated  the  DFIG  inside  wind  energy  conversion  system  using  the  classical  PI  and 
the  FOPI  controllers.  The  following  simulation  results  were  obtained  for  reactive  stator  power  reference 
Qs-ref  ~  0  anc*  wind  speed  profile  modeled  using  the  sum  of  harmonics  [21],  [22]: 

v  (t)  =  8  +  0.2  sin(0.10  t)  +  2  sin(0,35  t)  +  sin(1.235  t)  +  0.2  sin(3.5  t)  (34) 

Furthermore,  a  disturbance  of  500  A  is  applied  on  the  output  of  the  currents  plants  from  t  =  10  s 
on,  and  a  Gaussian  distributed  random  signal  measurement  noise  with  zero  mean  and  30  variance  value  is 
assumed  in  the  output  of  the  currents  plants  from  t  —  20  s  on. 

Figures  6(a),  (b),  (c)  and  (d)  present  the  evolution  with  time  of  the  wind  speed,  the  power 
coefficient,  the  wind  turbine  power  and  DFIG  slip,  respectively.  As  shown  in  Figure  6(b),  the  power 
coefficient  Cp  achieved  a  maximum  value  0.55  rapidly  and  remains  constant  over  wind  speed  evolution. 
However,  the  wind  turbine  power  and  the  slip  vary  according  to  the  wind  speed  profile. 

Electromagnetic  torque,  stator  active  and  reactive  powers  of  the  DFIG  are  shown  respectively  in 
Figure  7(a),  (b)  and  (c).  From  these  figures  it  is  clear  that  the  WECS  with  the  FOPI  controller  exhibits  better 
dynamic  performance  in  terms  of  the  set-point  reference,  the  disturbance  rejection  and  the  noise 
minimization,  than  the  WECS  with  classical  PI  controller. 
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Figures  6.  (a)  Wind  speed  profile,  (b)  power  coefficient,  (c)  wind  turbine  power  and,  (d)  DFIG  slip 


Time  (s) 


Figure  7.  (a)  Electromagnetic  torque,  (b)  stator  active  power,  (c)  stator  reactive  powers 


7.  CONCLUSION 

In  this  paper,  we  proposed  a  robust  control  of  a  variable  speed  DFIG-based  WECS  using  vector 
control  strategy  based  on  two  Fractional  Order  PI  Controllers,  to  prevent  system  deterioration  under  critical 
conditions,  namely:  DFIG  parameters  variation,  disturbance  and  measurement  noise.  The  two  Fractional 
Order  PI  controllers  were  designed  with  a  simple  analytical  method.  The  simulation  results  show  that  the 
WECS  with  the  FOPI  controllers  exhibit  better  dynamic  performance  even  in  critical  conditions.  This  is 
reflected  in  the  robustness  to  rotor  resistance  variation,  superior  disturbance-rejection  and  measurement  noise 
minimization. 
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